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Abstract. A stationary self-consistent outflow of a mag- 
netised relativistic plasma from a rotating object with an 
initially monopole-like magnetic field is investigated in the 
ideal MHD approximation under the condition <j/Uq > 1, 
where a is the ratio of the Poynting flux over the mass en- 
ergy flux at the equator and the surface of the star, with 
Uq = 7oWo/c and 70 the initial four- velocity and Lorentz 
factor of the plasma. The mechanism of the magnetocen- 
trifugal acceleration and self-collimation of the relativistic 
plasma is investigated. A jet-like relativistic flow along 
the axis of rotation is found in the steady-state solution 
under the condition cf/Uq > 1 with properties predicted 
analytically. The amount of the collimated matter in the 
jet is rather small in comparison to the total mass flux in 
the wind. An explanation for the weak self-collimation of 
relativistic winds is given. 
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1. Introduction 



Radio pulsars (Michel 1991), superluminal galactic objects 
(Mirabel & Rodri guez |1994| , |1998D and AGN s (Begelman, 
Blandford & Rees 1984, Pelletier et al. 1996| ) are observed 
to be associated with relativistic outflows. AGNs and su- 
perluminal galactic objects produce jets, while radio pul- 
sars demonstrate a surprisingly high efficiency of trans- 
forming the pulsar rotational energy into kinetic energy 



of the relativistic wind (Kennel & Coroniti 1984). Hence, 
the problem of acceleration and collimation of relativistic 
plasmas is among the most important in the physics of 
relativistic outflows. 

Nevertheless, radio pulsars, AGNs and superluminal 
galactic objects are objects of rather different nature and 
scales. It may well be that the mechanisms of acceleration 
of the outflows in these sources are not the same. How- 
ever, there seems to exist a mechanism which operates 
very well for all these objects, namely, the so-called mag- 
netocentrifugal acceleration where the acceleration of the 
plasma outflow occurs in the magnetic field of the rotat- 



ing object. Such conditions are realized in the magneto- 
spheres of radio pulsars, in the winds from accretion disks 
(Blandford & Payne 1982) and apparently in superluminal 
galactic objects (Mirabel & Rodriguez 1998). However, it 
has not been investigated whether the magnetocentrifugal 
mechanism could provide the acceleration of the relativis- 
tic plasma in these objects to the observed bulk motion 
velocities. 

The analysis of problem of the relativistic plasma ac- 
celeration and collimation in a rotating magnetosphere 



started with Michel (1969), where it was shown that the 



magnetocentrifugal acceleration of the plasma is not effec- 
tive for a prescribed monopole-like magnetic field. Subse- 
quently, the relativistic plasma outfiow has been investi- 
gated in several other studies and for di fferen t astrophysi- 
cal conditions (Okamoto 1978, Ardavan 1979, Beskin et al. 



1983, Camcnzind 1986). In particular, the ultrarelativistic 



plasma outflow from an accretion disk was investigated 
numerically in Camenzind ( 1987 ) while the structure of 
electromagnetically driven relativistic jets was discussed 
in Li ( 1993| ) under the assumption of jet confinement by an 
ambient medium. Collimation of stellar magnetospheres to 
jets in the force-free approximation was studied in Sulka- 
nen & Lovelace (1990|) and Fendt et al. (1995). Later, 



Fendt & Camenzind ( |1996 ) investigated the plasma dy- 
namics in the MHD approximation in magnetic configu- 
rations obtained via the force-free approximation. A class 
of self-similar solutions of the relativistic plasma outflow 
has been obtained in Contopoulos (1994). It was found 
that the efficiency of the plasma acceleration can be higher 
than in Michel's solution due to the effect of the magnetic 
collimation of the wind to the axis of rotation (Begelman 
& Li (|199^ ), Fendt & Camenzind |l996| , Takahashi & Shi- 
bata (1998)). Thus, the problem of the relativistic plasma 



acceleration in a rotating magnetic field appeared to be 
closely connected to the problem of plasma collimation. 

Strong evidence has been accumulated for a correlation 
between collimated astrophysical outflows (jets) and ac- 
cretion disks. Direct evidence for a disk-jet connection has 
been demonstrated by HST observations of the Herbig- 
Haro object HH30 (Burrows et al. 1996), although on the 
theoretical front this connection was discussed earlier (Ca- 
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Ferreira & Pelletier 1993, Falke & Bier 



1995). A disk-jet connection is also supported by 



observations of the active galaxy M87 (Ford et al. 1994) 
and superluminal galactic objects (Mirabel & Rodriguez 
199S). These observations seem to lead to the conclu- 



sion that jets require the existence of accretion disks for 



their formation (e.g. Livio 1999) and therefore accretion 
disks should be a necessary element of the collimation 
mechanism. However, we have by now reliable evidence 
that there is at least one source producing jets without 
the presence of an accretion disk. Namely, recent obser- 
vations of the Crab Nebula by the Chandra observatory 
leave no doubts t hat th e Crab pulsar indeed produces jets 
(Weiskopff et al. ^000 ). Therefore, it would be important 
to understand the mechanism of plasma collimation in all 
known jet sources, including the Crab pulsar. 

Such a mechanism is provided by magnetic collimation 
of outflows fro m acc retion disks a nd it w as first suggested 
by Lovelace (1976), Blandford ( 1976| ) and Bisnovatyi- 
Kogan & Ruzmaikin (1976). Subsequently, it was shown 
from an analysis of the full set of the MHD equations for 
nonrelativistic (Heayvaerts & Norman 1989) and relativis- 



tic outflows (Chiueh et al. |1991| , Bogovalov [1995| ) that col- 
limation is a rather general property of magnetised winds 
from a rotating object. A magnetised fast-mode supersonic 
wind inevitably will be partially cylindrically collimated 
along the axis of rotation at distances large compared to 
all characteristic dimensions of the source, such as the 
size of the accretion disk, provided that: the flow is non- 
dissipative, the magnetised central object is rotating, the 
total magnetic flux of one polarity reaching infinity (open 
field lines) is finite, the polytropic index of the plasma 
exceeds unity and the pressure of ambient matter is low 
enough not to be able to terminate the wind before the 
formation of the jet. This result is independent of the spe- 
cific properties of the central source. Collimation occurs 
spontaneously due to the internal Lorentz force arising in 
the magnetised plasma which is ejected from the rotating 
source and therefore can be regarded as a process of self- 
collimation (hereafter SC process). Note that the assump- 
tion of a polytropic gas with a constant polytropic index 7 
does not seem really restrictive, since several exact MHD 
solutions for collimated outflows have been constructed for 
the case of a variable index 7, Sauty & Tsinganos (1994), 
Vlahakis & Tsinganos (1998). 

Altho ugh t he SC process has been questioned (e.g., 
Okamoto 1999| ), a numerical investigation of this process 
and the characteristics of jets produced via this mecha- 
nism in the cold (Bogovalov & Tsinganos |1999D or hot 
(Tsinganos & Bogovalov pOOOl , Sauty et al (1999)) nonrel- 
ativistic winds leave no doubt that this mechanism works. 
This mechanism provides the formation of cylindrically 
collimated outflows at a large distance from the central 
source while collimation is accompanied by essentially 
more efficient acceleration of the plasma at the equato- 
rial region as compared with the original study of Michel 



( 1969| ). Evidence for the focusing of MHD outflows from 
accretion disks to the axis of rotation at distances of the 
order of the dimension of accretion disks has been ob- 
tained in other studies as well, such as in Ouyed & Pudritz 



(1997), Ustyugova et.al. (1999) and Krasnopolskii et. al. 
( |l999^ . 

In the analysis of Bogovalov & Tsinganos (1999), the 
relaxation method was used to obtain the solution of the 
steady state outflow problem, i.e., the time-dependent 
problem was solved in order to get a steady-state solu- 
tion. This method was first employed by Washimi & Shi- 



bata (1993) for outfiows from stellar objects. Later this 



method was applied in Bogovalov (1996) for the inves- 
tigation of the SC process in an initially monopole-like 
magnetic field. It seems that this numerical method al- 
lows us to solve accurately all problems connected with 
the c ritical surfaces (Bogovalov 1997b , Tsinganos et al. 



19961) 



In our previous study (Bogovalov 1997a), the relax- 



ation method has been applied to the investigation of the 
relativistic plasma outflow and it was found that the rela- 
tivistic plasma is very weakly collimated and accelerated 
even if the Foynting flux dominates the matter energy flux 
at the equator. However, an analytical study performed in 
this work has shown that an efficient SC of the relativistic 
plasma is still possible for very fast rotation of the cen- 
tral source. Recently it has been argued that a relativis- 
tic plasma is poorly accelerated even under this condition 
(Beskin et al. 199^ ). However, the SC process was totally 
neglected in this work. Therefore, the question concerning 
the efficiency of the acceleration and SC of the relativistic 
winds from fast rotators remains open. 

2. Basic assumptions 

Outflows of magnetised plasmas are described by the sys- 
tem of the familiar nonlinear MHD equations. It is not sur- 
prising that this system has many solutions corresponding 
to different physical conditions. Then, a reasonable choice 
of the appropriate approximation and a sound model is 
crucially important to arrive at some physically meaning- 
ful solution. 

To formulate the basic demands of the model let us 
consider the general properties of winds and jets. Accord- 
ing to observations, the length of the jets greatly exceeds 
all characteristic scales of the central source. Therefore, 
for an investigation of the SC process and the character- 
istics of the jets, the theory should allow us to obtain the 
solution on scales greatly exceeding all sizes of the central 
source. 

Observations of winds and jets from several classes of 



astrophysical objects such as our Sun (Parker 1963), or, 
YSO's (Burrows et al. 1996| ), pulsars (Kennel & Coroniti 



1984) and AGNs (Ferrari et al. 1996) show that they are 



super-magnetosonic at relatively large distances from the 
source. At the same time all winds must be subsonic near 
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the central source where the basic parameters of the out- 
flows are controlled. Therefore, the theory must describe 
flows of mixed type, containing the small scale regions of 
a subsonic flow and the large scale regions of the super 
magnetosonic flow. The scales of these regions may dif- 
fer by several orders of magnitude. These properties of 
the outflows dictate the choice of the approximations and 
mathematical models. 



2.1. Approximations 

The force-free approximation has been widely used for 
the analysis of axisymmetric relativistic outflows (Scharle- 
mann . 



Wagoner |1973| , Mestel & W ang |1979i Beskin et al. 
1983| , Lyubarskii |1990| ' Fendt et al. |1995D . For simplifica- 



tion reasons, the axisymmetric stellar magnetosphere can 
be assumed ini tially to have a dipolar magnetic field (Con- 
topoulos et al. 1999). However, this approximation is not 
appropriate for the investigation of the SC process for sev- 
eral reasons. For example, it describes only subsonic flows 
of plasma since the force-free transfield (Gr ad-Sh afranov) 
equation is elliptic everywhere (Fendt et al. 1995| ). There- 
fore, this approximation cannot be used for the investi- 
gation of plasma outflow in the supersonic regions. More- 
over, the coUimating and accelerating forces are equal to 
zero in this approximation. For these reasons, the force- 
free approximation gives ambiguous results. For example, 
the structure of the flow of the relativistic plasma at large 
distances from the source obtained in Sulkanen & Lovelace 



199C and in Fendt et al. 1995 differs from the structure 



obtained in Contopoulos et al. (1999). Although Sulkanen 
& Lovelace (1990) and Contopoulos ct al. (1999) used sim- 
ilar models, Contopoulos et al. (1999) obtained a radially 
outflowing plasma while Sulkanen & Lovelace (199C) ob- 



tained coUimated jets. Jets were obtained as well in Fendt 



et al. (1995). Actually all these works are mathematically 



correct. The difference between them is only in the bound- 
ary conditions at inflnity or at the outer boundary of the 
integration domain. These works clearly demonstrate that 
the flow of plasma at large distances from the source in the 
force-free approximation basically depends on the bound- 
ary conditions at infinity. But there is no such dependence 
in real supersonic flows. This is the most important quali- 
tative difference between the subsonic and super fast mag- 
netosonic flows. 

In the steady state subsonic flow the MHD signals 
propagate down and upstream of the flow. That is why the 
plasma flow is defined by the forces affecting the plasma 
at a given point, upstream and down stream of the fiow. 
This is not the case in the super fast magnetosonic fiow. 
No MHD signals can propagate upstream of the super- 
sonic flow. The supersonic fiow is defined only by the forces 
which affect the plasma upstream and at the given point 
and does not depend on the forces which affect the plasma 
downstream of the flow. For this reason the supersonic 
flow does not depend on the boundary conditions located 



downstream of the flow (actually the situation is more 
complicated; the flow should cross the fast mode critical 
surface, Bogovalov 1997h| ). The simplest approximation 
which allows us to describe the mixed type outflows is the 
ideal MHD approximation. This approximation is used in 
this paper. 

2.2. Model 

However, even the use of the MHD approximation does 
not always allow us to investigate the SC process in the 
supersonic region. There is a wide class of self-similar so- 



lutions for the nonrelatiyistic (Blandford & Payne 1982 



Vlahakis & Tsinganos 1999 e.g.) and relativistic plasma 
(Contopoulos 1994). All these solutions terminate not so 
far from the source. Independent of why this termination 
happens, the self-similar solutions cannot be used for the 
verification of the predictions of the SC theory since the 
central source has infinitely large size with infinite mag- 
netic flux in the upper hemisphere. At any distance from 
the source, the dimension of the source exceeds the dis- 
tance to the source. However, we are interested in the 
opposite case, when the distance is large compared with 
the dimension of the central source. Therefore, we demand 
that the appropriate model should have finite geometrical 
dimensions, finite magnetic flux of the open fleld lines of 
every polarity and initially uncoUimated plasma outflow. 
The last condition is important for us since our objective is 
to study the process which transforms the initially uncol- 
limated outflow into the outflow containing the jet. The 
word "initially" means that if one of the conditions for 
the SC process is not fulfilled (for example if the central 
source is not rotating or the magnetic field is absent) the 
wind expands radially at large distances from the source. 
Actually, a lot of models satisfy all these conditions. 

The SC is a general property of magnetised winds from 
rotating sources. It must take place in a wind from any 
source provided that the conditions specified by Heay- 
vaerts & Norman ( |1989| ), Chiueh et al. ( |199lD and Bo- 
govalov (1995) are fulfilled. Therefore "convenience" and 



"simplicity" become reasonable arguments for the choice 
of the model for the investigation of this process. The 
model of a star with a wind in an initially monopole- 
l ike (b elow m-1) magnetic field firstly proposed by Sakurai 
( 1985 ) is the most simple and convenient. This model al- 
lows us to avoid the mathematically complicated problem 
connected with the coexistence of open and closed field 
lines near the source which actually has no direct rela- 
tion to the SC process and usually is considered particular 
(Keppens & Goedbloed [T999|) . 



The model with the initially monopolc-like magnetic 
field describes more a wider class of outflows, including 
the time-dependent outflows from oblique rotators (Bo- 



govalov 1999) and outflows from objects with arbitrary 



distribution of the polarity o f the magnetic field at the 
base (Tsinganos & Bogovalov ^OOOl) . 
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In this paper we consider the outflow of an initially 
cold relativistic plasma and neglect only gravitation since 
we are interested in the effect of the electromagnetic forces 
only on the plasma. 

3. Basic equations and methods of solution. 

The plasma outflow from astrophysical objects consists of 
two zones where the properties of the flow differ qualita- 
tively. The nearest zone is located near the surface of the 
central source and includes all the critical surfaces of the 



steady state flow (Bogovalov 1997b ). The outer boundary 
of the nearest zone can be located rather arbitrary but 
always downstream of the fast mode critical surface. The 
flow of the plasma in the nearest zone is of the mixed 
type. It is subsonic near the source and is super fast mag- 
netosonic at the outer boundary. The far zone is located 
downstream of the nearest zone. The plasma flow in the far 
zone is super fast magnetosonic. Different methods should 
be applied for the solution of the problem in each of these 
zones. 

The solution of the problem of the steady state 
plasma outflow can be divided in two steps (Bogovalov & 
Tsinganos 1999). In the first step, the problem is solved in 
the nearest zone by the relaxation method and in the sec- 
ond, the steady state solution in the far zone is obtained by 
solving directly the transfield equation as a Cauchy prob- 
lem. The solution in the nearest zone is used to specify 
the initial values for the flow in the far zone. 

3.1. Solution in the nearest zone. 

The system of equations defining the dynamics of the rel- 
ativistic plasma in the ideal MHD approximation is as 
follows. 



pc(U + (v • V)U) = gE + ij X B, 

c 



E-l- -V X B = 0, 



^ ^ 47r . 1 . 

V X B = — j + -E, 

c c 



-B = -V X E, 

c 

p + (V • vp) = 0. 



(1) 

(2) 

(3) 

(4) 
(5) 



Here g = ( V • E) /An is the induced electric charge den- 
sity, j is the electric current density, E is the electric field, 
U = V7/c is the spatial component of the plasma four- 
velocity and p is the mass density of the plasma. 

The problem was solved in the cylindrical system of co- 
ordinates. Some details of the solution method are given 
in appendix A. The numerical simulation was performed 
only in the quarter of the total box of simulation due to the 
symmetry of the flow in relation to the equatorial plane 



and the axis of rotation. The boundaries of the nearest 
zone include the surface of the star, the axis of rotation, 
the equator and the outer boundaries, which are a cylin- 
der (right hand side of the box) and a plane perpendicular 
the axis of rotation at some distance from the source (up- 
per side of the box). The outer boundaries are located in 
the supersonic domain of the flow beyond the fast mode 
critical surface. This guarantees the independence of the 
solution in the nearest zone of the position and conditions 
of the outer boundaries. 

The number of independent boundary conditions on 
the surface of the star is equal to the number of MHD 
waves leaving the surface to satisfy the causality in the 
flow (Bogovalov 1997b ). They are as follows: 

1. A constant plasma density po- 

2. A constant Lorentz-factor 70 of the plasma in the coro- 
tating frame system. 

3. A constant normal component of the poloidal magnetic 
field Bq. 

4. The continuity of the tangential component of the 
poloidal electric field. 

5. A constant temperature of the plasma equal to zero. 

The boundary conditions at the axis of rotation cor- 
respond to the axial symmetry of the flow. The boundary 
conditions at the equator follow from the symmetry of 
the flow in relation to the equatorial plane. This means 
that p{z) = p{—z), Vj,{z) = Vj,{~z), Vz{z) — -~Vz{—z), 
Br{z) = Br{-z), Bz{z) = -Bz{~z), v^{z) = v^{-z) and 
B^{z) = B^(-z). 

At the outer boundaries only one side internal deriva- 
tives were used for calculation of the derivatives in the 
finite differences, since the solution does not depend on 
the conditions down the flow from these boundaries. 

3.2. Solution in the far zone 

The flow in the far zone is fast mode supersonic. The sys- 
tem of equations describing the steady state flow is hyper- 
bolic. Therefore, a Cauchy problem can be formulated for 
this flow. 

It is convenient to introduce the flux function ijj as 
follows 



Bp — 



(6) 



Where Bp is the poloidal magnetic field of the axisymmet- 
ric outflow and ip is the unit azimuthal vector. 

The stationary MHD equations admit well known four 
integrals (Tsinganos 1982| ). They are: 

(a) The ratio of the poloidal magnetic and mass fluxes, 
(/3) The total angular momentum per unit mass L{ip), 



rU^ 



FrB^ = L{i,) . 



(8) 
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(7) The corotation frequency ^(ip) in the frozen-in con- 
dition 



U^Bp - UpB^ ^ rjBpn{'il;)/c. 



(9) 



{S) The total energy W{ip) per particle in the equation 
for total energy conservation, 



(10) 



This system of equations should be supplemented by a rel- 
ativistic relationship between the components of the four- 
velocity 



7 



1 + C/2 + (7^. 



(11) 



Momentum balance across the poloidal field lines is 
expressed by the transfield equation. For analysing the 
behavior of the plasma at large distances, it is convenient 
to deal with this equation in an orthogonal curvilinear 
coordinate system ('0,??) formed by the poloidal magnetic 
field lines and by the lines of the poloidal electric field. 
ip varies with the motion across the magnetic field lines, 
while r] varies with the motion along the magnetic field 
lines. A geometrical interval in these coordinates can be 
expressed as 



(12) 



where gip,gri are the corresponding line elements (compo- 
nents of the metric tensor). 

According to Landau & Lifshitz (1975) the equation 
T^.^, = 0, where T*-' is the energy-momentum tensor and 
"; k" denotes covariant differentiation, will have the fol- 
lowing form in these coordinates. 



_a_ 

dip 





— 1 9r 


8-K 


r dip 



Bt-E^ 



1 dgr, 
g-q dip 



UpVpCp - 



47r 







(13) 



with the electric field E — —B„. 

c t> 

The unknown variables here are zirj^tp) and r{rj^'ip). 
The metric coefficient grj can be obtained from the trans- 
field equation (|l^). 



V 

g^ = exp{J G{?], ip)dip) 



(14) 



where 



G(77,V') 



d 


B^-E^ 


_ 1 dr 


dtp 




r dip 



U^pV^cp - 



B^ -E^ 



UpVpCp - 



4tt 



(15) 



The lower limit of the integration in (|lj) is chosen to be 
such that the coordinate rj is uniquely defined. In this 
way r] coincides with the coordinate z where the surface 
of constant 77 crosses the axis of rotation. 



The metric coefficient g^ can be obtained from defini- 
tion (|^) in terms of the magnitude of the poloidal magnetic 
field as follows 

Sip ^ -ZTT ■ (16) 



rB 



p 



The orthogonality condition 

^rjl^Tp Z^Z^ — , 

and the relationships 



2 _ 2 I 2 

9ri ^ "^fj ' ■^rji 



9% 



2,2 

Tih + '^ib 



give us that 

^ip9v 

9ip 



ripgr, 
9ip 



(17) 
(18) 

(19) 



with 5^ calculated by the expression ([l4|). Here = 
dr/drj, = dz/drj, = dr/dijj, z^ — dz/dtp. For the 
numerical solution of the system of equations (^9|) the two 
step Lax-Wendroff method on the lattice with a dimension 
equal to 1000 is used. 

Equations (|l^) should be supplemented by appropri- 
ate boundary conditions and initial values on some ini- 
tial surface of constant rj located in the nearest zone, but 
down the flow from all the critical surfaces. The form of 
the initial surface of constant rj was obtained numerically 
through integration of the following equations. 



dr 

dip 



B, 



rBf 



dz 
9^ 



Br 



rBj 



(20) 



The values B^, Br and integrals 
W{il}),L{ip)^^{ip)^F{4>) were taken from the solution of 
the problem in the nearest zone. 

Boundary conditions on the axis of rotation and the 
equatorial plane are the same as the conditions in the 
nearest zone. No conditions at infinity were specified. 

4. Parameterisation of the solution 

A stationary solution of the relativistic plasma flow from 
the star with an initially m-1 magnetic field is d efined 
by three dimensionless parameters (Bogovalov 1999| ) . One 
of them is the ratio of the radius of the star Rstar to 
the radius of the initial fast mode magnetosonic surface 
Rf = ^ {BaRstarY /^t^PoVqcUo, whcrc C/q = «o7o/c. The 
word "initial" here refers to the radius of the fast mode 
surface when the star is not rotating. But this dependence 
can be neglected under the condition Rstar Rf and 
Rstar ^ Rii where Ri — c/fl is the radius of the fight 
cylinder, and f2 is the angular velocity of the uniformly 
rotating star. Therefore, the solution crucially depends 
only on two dimensionless parameters which can be cho- 
sen rather arbitrary. It is convenient to choose one of them 
as follows 



47rpo7oWo 



^ Rstar^ ^2 



(21) 
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(7 is the squared ratio of the initial radius of the fast mode 
surface to the light cylinder radius. According to the 
radius of the light cylinder in units of i?/ is 



3l 



(22) 



The parameter a is closely related to the so-called mag- 
netisation parameter usually used in the physics of the 
relativistic winds from radio pulsars (Arons 1996| ). The 
magnetisation parameter is the ratio of the Poynting flux 
density to the flux density of the mass energy (including 
rest mass energy). The parameter a used here is equal to 
the ratio of the Poynting flux to the mass energy flux at the 
equator of the rotator with the initially m-1 magnetic field 
provided that c/C/g << 1 (Bogovalov 1999). In our pre- 



vious work (Bogovalov 1997a), we used another definition 
of the a parameter introduced initially by Michel ( 1969| ). 
This parameter is denoted by (Tj\/. There is a simple re- 
lationship between them am ~ (tUq. It follows from this 
relationship that cta/ is equal to the Poynting flux density 
per particle on the equator provided that <j/Uq << 1. 
It is convenient to take the second parameter as follows 



Rf 



RiUo 



(23) 



The physical sense of this parameter becomes clear if we 
pay attention to the fact that in the model of Michel 
( |1969D the terminating four-velocity of the plasma on the 

1 /3 

equator is Uoc = {ctat ) and a can be presented as follows 



(24) 



According to Eq. (|24|), the parameter a shows how effec- 
tively the magnetic forces arising in the rotating magnetic 
field accelerate the plasma. If a < 1 acceleration is not 
important compared to the initial velocity of the plasma. 
It is natural to call the rotators with a < 1 as slow rota- 
tors. Correspondingly, the rotators with a > 1 should be 
fast rotators which could in principle (but not necessar- 
ily) accelerate the plasma. The parameter a introduced 
here is the relativistic generalization of the a parameter 
introduced earlier by Bogovalov & Tsinganos (1999) for 
the nonrelativistic plasma. Numerical simulations of the 
nonrelativistic plasma outflow show that this parameter 
indeed divides rotators into two groups with different effi- 
ciency of acceleration and collimation of the plasma. Be- 
low we will see that this division remains valid also in the 
relativistic limit. 

Representation of the solution in the parameters a, cr is 
convenient since the steady state outflow in the initially m- 
1 magnetic field has a simple analytical solution in the limit 
a independent of the value of a (Bogovalov 1997a) 



and the transition to the nonrelativistic limit corresponds 
simply to the limit cr ^ and uq/c <C 1 . 

For comparison with observation, Eq. (|2l| ) should be 
expressed in variables general for the m-1 model and real 



astrophysical objects. An astrophysical object can be con- 
sidered as equivalent to the object with the m-1 magnetic 
field if they both produce winds with the same magnetic 
flux, the same mass flux, the same initial Lorentz factor 
and finally both rotate with the same angular velocity. 
The magnetisation parameter in these variables takes the 
form 



(25) 



where ipt is the open field lines flux of one polarity in 
the wind, N is the rate of ejection of particles with 
mass m. For radio pulsars, ■0* = T^BoRp where Rp = 
Rstar \/ Rstar^/c is the radius of the polar ca p on t he sur- 



face of the pulsar ( Ruderman & Sutherland |1975[ ) . 

The position of the radio pulsars in the parameter 



space a, cr was discussed in detail by Bogovalov (1999). 
All radio pulsars eject Poynting-dominated winds (cr > 1) 
but basically all of them are slow rotators (a < 1). One 
of the youngest and fastest rotating pulsars with a strong 
magnetic fleld on the surface is the pulsar in the Crab 
Nebula. This pulsar rotates with Q. 200 and ejects 
pairs with the rate iV « 4 • 10'^® per s econd and average 
70 ~ 200 ( Daugherty & Harding |1982|) . At Bo « 4 • lO^^ 
G the magnetisation parameter of the wind cr w 1.3 • 10"* 
and a « 0.57 The light cylinder radius of the Crab pul- 
sar is located at the distance Rl = 1.5 • 10* cm from the 
axis of rotation. The radius of the fast magnetosonic sur- 
face Rf exceeds the light cylinder radius by a factor of 
114. Marginally, the Crab pulsar could be the fast rotator 
provided that 70 < 200. The uncertainties in estimates 
of this parameter do not rule out this possibility. But all 
other, more slowly rotating pulsars are certainly the slow 
rotators. 

In AGNs and superluminal sources 70 is deflned from 
observations of jets and lies in the range 1 — 10. If the 
collimation of the plasma is really due to the SC process 
we should expect that a > 1 and a > Jq. So, the winds 
from AGNs and superluminal sources should be mainly 
Poynting flux dominated. 

5. Basic results 

In Fig. we show the result of the numerical simulation 
of the plasma outflow in the nearest zone of the star, for 
a time interval t = 2.62Rf/vQ and for the set of the values 
of the two key parameters, a = 1.0 and a = 51. This time 
interval is larger by a factor of 2.62 compared to the typ- 
ical time the plasma and all MHD perturbations take to 
travel from the surface of the star to the fast mode surface. 
Note that during this time interval, a steady state solu- 
tion is reached, [] since all perturbations which could be 

This solution is stable. In the previous work (Bogovalov 



1997a) the time was expressed in the units 1/n. Therefore, 



the ratio of the simulation time in the present work over the 
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0=5 1 a=l.O Yq=7 dim=200x200 t=2.62 




Fig. 1. The steady-state relativistic plasma flow in the 
nearest zone of the star for the parameters ct, a, initial 
Lorentz-factor 70, dimension of the lattice and duration 
of the simulation in units Rf/vo indicated at the top of 
the figure. The cylindrical coordinates z (axis of rotation) 
and r (equator) are given in units of Rf. The star is lo- 
cated in the lower left-hand corner of the figure. Solid 
lines radiating out from the star indicate the lines of the 
poloidal magnetic field while dashed lines represent the 
lines of the poloidal electric currents generated at the ro- 
tation. Thick (thin) solid lines indicate the fast (Alfvenic) 
mode surface, while the vertical dot-dashed line represents 
the light cylinder. 



excited at the start of the simulation have left the box of 
simulation. The initial fast mode surface is circular with a 
radius equal to 1 in our units. Note that the final fast mode 
surface shown in Fig. 1 has a small amplitude kink struc- 
ture at a polar angle around 45 degrees. Such a structure 
is formed in all simulations performed with the present 
method and in the present geometry (Washimi & Shibata 
Tsinganos 19991 ). Basically it is due 



1993, Bogovalov 
to the shape of the surface of the star which is actually 
approximated by a sequence of step functions in the cylin- 
drical system of coordinates. 



time in the previous work is 



present / ipreviou. 



:){Rfn/vo) 



{tpresent/tprevious)\i^ provided that VQ ~ c. Here t 



present 



and 



tprevious are the dimensionless times in the present and pre- 
vious works. In the present work the solution is stable up to 
the time tpresent = 2.62 at a = 1 and a = 51. In the previous 
work the solution for a — 0.0014 and a — 3.3 was destroyed 



after t„ 



= 2.02. Thus, in the present work much stronger 



Poynting-dominated wind was simulated, for 9.26 times longer 
than in the previous one. The simulation shown in the Fig. 4 
lasted even longer. This confirms the stability of the solution. 
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Fig. 2. The relativistic plasma flow in the far zone for 
a — 1.0 and ct = 51. The coordinates are expressed in 
units of Rf. 

The SC theory predicts that a wind of relativistic 
plasma consists of a cylindrically collimated core sur- 
rounded by a radially expanding flow (Heyvaerts & Nor- 
man |l989[ Chiueh et al. [1991| and Bogovalov p^95| ) . An in- 
spection of Fig. 1^ shows that the collimation of the plasma 
along the axis of rotation is practically invisible in the 
nearest zone. The results of the same simulation in the far 
zone are presented in Fig. At a flrst glance, in the far 
zone there is no cylindrical collimation. However, a closer 
examination of the results shows that a cylindrically col- 
limated core is actually formed in the wind. 

Fig. ^ gives the distribution of the poloidal magnetic 
field across the axis of rotation for two values of z. This 
distribution differs from that expected for a m-1 wind. The 
m-1 poloidal magnetic fleld is described by the equation 
Bp = Bo/{r'^ + z"^). And, at a fixed z the poloidal field nor- 
malized to its value at the axis of rotation (Baxis) depends 
on r as Bp/Baxis = 1/(1 + {r/zf). Then, at z = 1121% 
the poloidal magnetic field remains practically constant 
atr<5^z = 1121. This dependence is represented 
by the horizontal line Bp/Baxis = 1 in Fig. ^. However, 
our numerical simulation gives a totally different picture. 
Namely, the poloidal magnetic field has clearly a maxi- 
mum at the rotation axis and then drops down with the 
distance r from the axis, for example at r ~ i?/ to half its 
value from the axis. It is important that this distribution 
is practically independent of z. This result implies that the 
wind near the axis of rotation has achieved an asymptotic 
state and the flow is collimated exactly along the axis of 
rotation here. 

The m-1 model allows us to perform a quantitative 
comparison of the SC theory and our numerical experi- 
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ment. The structure of the cylindrically coUimated flow 
at large distances from the source is especially simple pro- 
vided that the central source rotates uniformly while the 
plasma velocity and the ratio of the poloidal magnetic and 
mass fluxes (function cF in Eq. (^) are constant across 
the jet. In this case the dependence of Bp on r at the 
fixed z for a cold plasma is given by the simple analytical 



expression (Bogovalov 1995) 



Bp 

B'n 



{r/R,r 



(26) 



where Rj = "fjVj/il = RfUj/Uoa is the characteristic ra- 
dius of the coUimated fiow of the cold plasma, while "fj and 
Vj are the Lorentz factor and the velocity of the plasma in 
this fiow, respectively. Eq. (26) is a direct consequence of 
the SC theory and is valid a,i r ^ Ri. The dependence of 
the normalised plasma density on r is identical to the dis- 
tribution ( p6| ) since F and the plasma velocity are assumed 
to be constant across the cylindrically coUimated flow. In 
this case the plasma density is simply proportional to the 
poloidal magnetic fleld, since p = Bp/4TTVpcF{ip) (see Eq. 
(0))- 

All these conditions are indeed fulflUed in our case. The 
outflow speed at the vicinity of the axis of rotation is close 
to the speed of light ( Vp = 0.87c at 70 = 2 and Vp = 0.98c 
at 7o = 7) and is independent of r. The function F is 
assumed independent of ip for all fleld lines and the star 
rotates uniformly. A comparison of the calculated depen- 
dence of Bp on r with the analytical prediction, Eq. (p6|), 
shows that the SC theory correctly predicts the asymp- 
totic properties of the wind structure. 

The cylindrically coUimated flow in Fig. (^) is not vis- 
ible since its transversal scale is extremely small. The flux 
of the poloidal magnetic fleld contained in this flow is small 
as well. The dependence of the magnetic flux on r is shown 
in Fig. (^) by the dashed-dotted line. For plotting conve- 
nience, this flux is multiplied by 100 and is also normalised 
to the total magnetic flux in the upper hemisphere, ipeq, 
which coincides with the magnetic flux at the equator. 
The mass flux normalised to the total mass flux in the up- 
per hemisphere is presented by the same curve since the 
mass flux is proportional to the poloidal magnetic fleld in 
our case. The fraction of the magnetic and mass flux in 
the cylindrically coUimated flow is of the order 0.1%. The 
rest of the magnetic and mass flux expands radially, in 
accordance with the SC theory. 

The same physical picture emerges for a larger value 
of a. In Fig. (^ we show in the nearest zone the results 
of the simulation of the fiow of plasma for a — 4. The 
increase a does not change significantly the structure of 
the flow, although collimation now becomes visible even in 
the nearest zone. The same flow in the far zone is shown in 
Fig. (^. It may be seen from this flgure that a very narrow 
coUimated beam is formed in the wind. The structure of 
the cylindrically coUimated flow is shown in Fig. (H) for 
two values of z. 




Fig. 3. The distribution of the poloidal magnetic fleld 
across the jet a.t z = 1121Rf (solid line) and z = 140i?/ 
(dots), compared with the prediction of analytical theory 
(dashed line) . The poloidal magnetic fleld is normalized to 
its value at the axis, r = 0. The dashed-dotted line gives 
the distribution of the ratio of the magnetic flux to the 
total magnetic flux in the upper hemisphere multiplied by 
100. The distance to the axis of rotation r is expressed in 
units oi Rf such that the radius of the jet is 1/a = 1. 

The efflciency of the acceleration of the plasma by a 
rotating object is important in astrophysical applications. 
This acceleration in the model under consideration is ex- 
pected to be most efficient near the equatorial plane since 
at the axis of rotation there is no accelerating force at all. 
The accelerating force in the poloidal direction is propor- 
tional to the product ji_B^, where j_l is the component of 
the electric current perpendicular to a field line. The ac- 
celerating force in the azimuthal direction is proportional 
to the product j±Bp. Thus, all accelerating forces are pro- 
portional to j_L . This implies that the acceleration occurs 
only if the electric currents cross the poloidal field lines. 
The electric currents shown by the dashed line in Figs. (|l]) 
and (||) indeed cross the poloidal field lines. This is very 
well seen especially in Fig. (^) near the equator. 

The dependence of the four-velocity of the plasma on 
the distance from the centre of the source is shown in 
Fig. (0). The distance is given in a logarithmic scale. The 
energy of the particles increases slower than logarithmi- 
cally and reaches a constant value. It follows from en- 
ergy conservation (|l^) that the maximum possible termi- 
nal four-velocity Uoo at the equator is Uoc = Uq(vooIvq)<j, 
since W at the equator can be estimated as 7o<t, pro- 
vided that (T ^ 1 and the toroidal magnetic field is 
estimated as B^ k, —rilBp/vo (Bogovalov 1997a). All 



Poynting fiux should be transformed into kinetic energy 
of the plasma in this case. But this value of the veloc- 
ity is never achieved. Michel ( |1969D calculated for Uoo the 
value Umichei = (ff/o)^^^ for the plasma fiow in a pre- 
scribed and fixed m-1 poloidal magnetic field. Therefore, 
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dim=300x300 t=3.34 




Fig. 4. Steady-state relativistic plasma flow in the near- 
est zone of the star and for the parameters a, a, initial 
Lorentz-factor 70, dimension of the lattice and the du- 
ration of the simulation in units Rf/vo indicated at the 
top of the figure. Coordinates z (axis of rotation) and r 
(equator) are given in units of Rf. The star is located in 
the lower left-hand corner of the figure. Solid lines radi- 
ating out from the star indicate the lines of the poloidal 
magnetic field while dashed lines represent the lines of 
the poloidal electric currents generated at the rotation. 
Thick (thin) solid lines indicate the fast (Alfvenic) mode 
surface, while the vertical dot-dashed line represents the 
light cylinder. 
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Fig. 5. Relativistic outflow in the far zone of the star for 
a = 4 and a = 48. The coordinates are expressed in units 
ofi?/. 
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Fig. 6. The dependence of the poloidal magnetic field 
across the jet at z = 10350i?/ (solid line) and a,t z = 
lOOOi?/ (dotted line) compared with the prediction of an- 
alytical theory (dashed line). The poloidal magnetic field 
is normalized to its value at the axis, r = 0. The dashed- 
dotted line gives the distribution of the ratio of the mag- 
netic flux to the total magnetic flux in the upper hemi- 
sphere multiplied by 100. The distance to the axis of ro- 
tation r is expressed in units oi Rf such that the radius 
of the jet is 1/a = 0.25. 



Umichei/Uo = c?l^ (see Eq. U). It is reasonable to esti- 
mate the efficiency of the plasma acceleration as the ratio 
C = U 00 1 U^a which is approximately the ratio of the part 
of the Poynting fiux transformed into the kinetic energy 
of the plasma (this is valid, provided that Uao ^ Uq) 
to the initial value of the Poynting flux at the equator 
in the relativistic limit. For the Michel solution ( 1969[ ) 
^ = c?!'^ ja. The same value is given by Beskin et al. 
(1998) for the analytical solution of the same problem un- 
der the assumption that the poloidal magnetic field goes 
to a monopole-like field at a ^ cxd. It follows from Fig. 
(R) that Uoc/Uq increases with a essentially faster than 



/3 



a"''^. For the fiow with a = 4 Michel's estimate gives 
Up.Qo/Uo = 2.5, while the numerical solution gives more 
than 8. Correspondingly, the efficiency of the acceleration 
remarkably exceeds the efficiency of the acceleration given 
by Michel and achieves 16.7% for the flow at a = 4 and 
cr = 48. It follows from Fig. |7| that the remarkable part 
of the acceleration occurs at large distances from the star. 
It happens at the collimation of the magnetic flux to the 
axis of rotation due to disbalance arising between the ten- 
sion and pressure gradient of the toroidal magnetic field 
(Begelman, & Li |1994 Fendt & Camenzind |1996| ). This 
mechanism was not taken into account in the analysis of 
Michel ( |1969| ) and Beskin et al. ( |1998D since they assumed 
a fixed m-1 magnetic field. 

It follows from Fig. ^that the largest energy is achieved 
not at the equator but at r = lOOOi?/ for a = 4. To 
understand why this happens it is convenient to look on 
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Fig. 7. Dependence of Up along a field line on the dis- 
tance to the centre of the source expressed in Rf. Curve 
1 is for a = 1.0, a ~ 51 for the equatorial plane. Curves 
under number 2 are for the equatorial plane (solid line), 
for ip/tpeq = 0,7 (dashed line), for ^p/tpeq = 0,3 (dashed- 
dotted line) and for ip/i^eq — 0,03 (dotted line) for the 
parameters a — A, a — A8. 

the dependence of Up/Uo on ip across field lines. This de- 
pendence is presented in Fig. (||) by two upper curves 
1 ( for z = lOOOi?/) and 2 (for z = 10350i?/). These 
curves have small amplitude irregular structure which is 
purely of numerical origin. It arises at large distances 
where the program deals with large numbers. The max- 
imum of the acceleration is achieved not at the equator, 
but at 'ip/ipeq ^ 0.3. This is explained by the fact that the 
asymptotical regime predicted by the SC theory is actually 
not achieved at these z for ipli'eq > 0.05. It appears that 
the larger ip the larger the distance which is necessary to 
reach the asymptotic regime of the outflow. This property 
of the SC process has already been found in the nonrela- 
tivistic outflow simulation (Bogovalov & Tsinganos p^99| ). 

The physics of this phenomena can be understood as 
follows. It can be shown from simple geometrical consid- 
eration that 



9v 9ip 



1 



rRrB. 



(27) 



p 



where Rc is the curvature radius of the poloidal magnetic 
field lines, with R^. positive if the centre of curvature is in 
the domain between the line and the axis of rotation and 
negative in the opposite case. With this expression, Eq. 
( p^ ) takes the form 



_a_ 

dip 



Sttt^ dip 



[r^{Bl-E^)\- 



1 

47rr 



dr 
dip 



UpF{iP) 



(28) 



[U^-F(iP){l-(rn/cf)B^] 
i-KrRaF(i>) 



0. 



Let us consider the flow in the region of radial ex- 
pansion, as in the case under consideration, where the 
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Fig. 8. Dependence of the four-velocity of the plasma 
Up/Uo (upper curves) and Heyvaerts & Norman integral 
Ihn (lower curves) on ip for z = lOOOi?/ (curve 1) and 
z = 10350i?/ (curve 2). The parameters of the outflow are 
a = 4, cr = 48 , 70 = 2. 

SC is small. For simplicity, we neglect the acceleration 
of the plasma. The first term in this equation decreases 
with r as r^*. The second one decreases as r"^, since 
E ^ rBp (Bogovalov 1997a). The third term is pro- 



portional to the product rU^Bp. This term decreases as 
r~^, since rBp ^ 1/r and U^p ^ 1/r due to angular mo- 
mentum conservation. Thus, only the second and fourth 
terms survive in Eq. (^8|) at these large distances. There- 
fore, this equation takes the form 



_l d_ 

Sttx"^ dtp ' 



{x\Bl-E^)) 



{Up + Fiib)x^Bp 
ATTxRcFi^P) 



0, 



(29) 



where x = r^l/c. From the frozen- in condition (Q) we 
may estimate the toroidal magnetic field at large distances 
when ^ Q, where is the azimuthal velocity. 



Bin — . 

V 



(30) 



In the asymptotic regime plasma moves along straight 
lines, when i?c — *■ cx). The structure of the flow in this 
regime is defined by the trivial equation 



d_ 

dip 



{x^Bl ^ E^)) ^ 0. 



(31) 



Together with Eq. ( pOj ) this equation gives that the rel- 
ativistic generalization of the Heyvaerts & Norman in- 
tegral Ihn = \xB^\/'-f (Chiueh et al. 1991) should be 
constant in the region of radial expansion. This integral, 
shown in Fig. (|^) by the two lower curves for essentially 
different distances, is actually constant in the small region 
0.01 < ip/tpeq < 0.05 (Heyvaerts & Norman plateau in the 
figure). But at tp/ipeq > 0.05, Ihn is not constant. This 
contradicts the expected behavior oi Imh, but, as we show 
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below, this implies that the asymptotic regime is still not 
reached in this region of ip. 

The characteristic scale of collimation Dcoii which de- 
fines the distance where the flow reaches the asymptotic 
regime can be estimated as follows. dO/dR = l/Rc, where 
9 is the polar angle of the plasma velocity and R is the dis- 
tance to the center of the source. In the radially expanding 
wind dip = BpT^ sm{6)dd. Taking into account also that 



at large distances — —{rfl/vp)Bp (Bogovalov 1997a), 
the rate of turning of the trajectory of the plasma follows 
from Eq. (E9|) and is defined by the equation 



l-l 



< 



sin 20 



(32) 



Integration of this equation gives the estimate of the dis- 
tance at which the line of the flow turns through the angle 
Ad. The scale of the collimation of the flow line DcoU is de- 
fined by the distance at which the flow line turns through 
the angle AO ~ 9. The integration of equation can be 
performed from the fast mode surface since this consider- 
ation is valid only in the supersonic region. The scale of 
the collimation is estimated as 



DcoU > Rf exp 



sm 



i29 



(33) 



According to this equation there are two regions of the 
flow with the essentially different dependence of Dcoii on 
the parameters of the flow. The smallest DcoU takes place 
for small angles limited by the value 9^ < 1/a. For these 
angles the scale of the collimation is (see also Bogovalov 
& Tsinganos |l99g|) 



DcoU > Rfcxp 



1 

2^' 



(34) 



Similar to the nonrelativistic outflow the collimation of 
the plasma close to the rotational axis occurs near the 
fast mode surface at a > 1. At the initially m-1 outflow 
the magnetic flux depends on the angle 9 as "tp — ''jjeqi^ — 
cos 9). Therefore, the part of the magnetic flux where the 
collimation occurs as in the nonrelativistic limit is 



1p 1 
— < — . 

- 2a 



(35) 



The collimation scale of the other part of the flow ( at low 
latitudes) is defined by the equation 



9 sin 9U^ 

DcoU > -R/ exp — ^. 

2 cos 9 



(36) 



and at the angles 9^ ^ '^/Uq the collimation scale becomes 
large compared with Rf independent of a. It follows from 
this result that the remarkable bend of the flow lines at 
low latitudes is possible only for the plasma with relatively 
low Uq. The collimation scale increases exponentially with 
the growth of the Lorentz factor of the plasma. These 
properties of the flow explain the numerical results. 



In the simulation we obtain such a behavior of the 
flow, with three regimes: a core where the current increases 
with ip, the Heyvaerts & Norman regime and the regime 
where the enclosed current increases (Fig. ^). Thus, an 
extremely narrow cylindrical core is formed which contains 
the magnetic flux ip < {l/2a)'tpeq- This core is formed at 
a distance comparable to Rf for a > 1, similarly to the 
nonrelativistic limit. For larger -0 the collimation scale is 
defined by Eq. ( |36| ) . For 70 = 7 the characteristic scale of 
coUimation is of the order of Rf ■ iQiO-^^*'^'!^. It imphes 
that the smaller 9, the earlier the asymptotic regime will 
be achieved. This explains why the Heyvaerts & Norman 
plateau is formed only at the relatively small 9 (small ip). 
At larger 9 the scale of collimation exponentially increases 
and finally diverges at 9 = tt/2. Thus, the asymptotic 
regime cannot be achieved at low latitudes. At the final 
distance from the source there should exist a region at low 
latitudes (large ip) where Inh is not constant. It is this 
behavior which is found in our calculations. 

For a = 4 the collimation is more efficient and the 
effect of collimation appears already in the nearest zone 
where the plasma is not yet accelerated. At the estimation 
of the collimation scale it is necessary to take into account 
acceleration of the plasma to a Lorentz factor greater than 
10 near the equator. This could be done by integration of 
Eq. ( p2| ) from a radius Rm > Rf where the energy of the 
plasma practically reaches the terminal value. If we take 
into account the acceleration, the asymptotic regime (in 
which Inh does not depend on -0 in a radially expanding 
wind) will be achieved at a distance of the order DcoU > 
^^-[^Q22 etane^ It follows from this expression that at a = 4 
the asymptotic regime at i? « 10^-i?/ will be achieved only 
for ip/ipeq < 0.1. This value agrees well with the numerical 
results presented in Fig. ^ in the limits of uncertainties. 

Comparison of the prese nt wo rk with our previous 
work (Bogovalov & Tsinganos 1999 ) shows that the rela- 
tivistic plasma is collimated less efficiently than the non- 
relativistic one. In particular, the fraction of the colli- 
mated mass flux in the nonrelativistic flow was close to 
1%, while in the relativistic limit it is of the order 0.1%. 
This difference can be understood from simple qualita- 
tive consideration. The dynamics of the cold plasma is 
basically controlled by two forces. One of them is the vol- 
umetric Lorenz force, (j x B)/c and the other is the vol- 
umetric Coulomb force, qE (see Eq.(l)). The collimation 
of the plasma to the axis of rotation is defined by the 
components of these forces perpendicular to the stream 
lines. The perpendicular component of the Lorenz force is 
Fl± = ^j^^Bip, where j|| is the component of the poloidal 
electric current parallel to the poloidal magnetic field and 
B^ is the toroidal magnetic field. The perpendicular com- 
ponent of the Coulomb force is Fc± = q ■ E± . These forces 
are shown schematically in Fig. (^. It is important to note 
that these two forces are directed along opposite direc- 
tions. The Lorenz force collimates the plasma while the 
Coulomb force decoUimates it. In the nonrelativistic limit 
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the Coulomb force is neglected, since it is negligible in 
comparison to the Lorenz force. The ratio of the Coulomb 
and Lorenz forces is of the order (w/c)^ in this Hmit. In 
the relativistic limit the Coulomb force becomes practi- 
cally equal to the Lorenz force and it strongly decreases 
the SC of the plasma. That is why the relativistic gen- 
eralization of the parameter a introduced earlier for the 
nonrelativistic winds (Bogovalov & Tsinganos 1999) has 
additionally the Lorentz factor in the denominator (see 
Eq. (H)). 

Although the Lorentz factor in the denominator of 
Eq. decreases remarkably the value of a, it was ex- 
pected that the collimation of the relativistic plasma will 
be strong enough at the condition a > 1. In the nonrela- 
tivistic case the collimation becomes strong already in the 
subsonic region under this condition. The most surpris- 
ing result is that the relativistic plasma is practically not 
coUimated, even under this condition. The increase of the 
collimating force does not provide more efficient collima- 
tion of the relativistic plasma! This happens since another 
relativistic effect comes into play. 

The curvature Kc = 1/^^c of the stream lines of the 
nonrelativistic plasma is defined by the simple equation 

F± 

Kc = 2 , (37) 

mnv 

where m is mass of particles, F± is the sum of all forces 
affecting the plasma perpendicular to the stream line. It 
is seen from this equation that an increase of the force 
F± results in an increase of Kc, such that collimation be- 
comes stronger. In other words, as it is evident from Eq. 
(^), Kc is proportional to the collimating force and in- 
versly proportional to the mass density of the outfiowing 
plasma. The situation changes essentially in the case of 
the relativistic, Poynting dominated plasma {a > 1) where 
the energy flux is dominated by the electromagnetic field. 
Then, according to the familiar relationship between en- 
ergy and mass {E = mc^), the flux of the energy is equiva- 
lent to some inertial mass flux. It follows from Eq.(p9|) that 
eq. ( ^7| ) for the relativistic Poynting-dominated plasma is 
modifled as follows 
F± 

Kc = ^2 — . (38) 

nv^^(m + -. — ^) 

This equation demonstrates that the effective mass of the 
Poynting- dominated plasma is dominated by the toroidal 
magnetic field. It is important that now the effective mass 
increases by a term proportional to S^. The collimating 
force is also proportional to the same term (see eq. (p9|)). 
Therefore, the increase of the collimating force is accom- 
panied by a corresponding increase of the effective mass 
density of the plasma. This in turn results in the fact that 
the curvature Kc of the stream lines does not increase with 
the increase of the collimating force. This is the reason 



6. Conclusion 

This work shows that the self-coUimation (SC) process oc- 
curs in relativistic winds as predicted by the SC theory. 
The wind consists of a radially expanding part which sur- 
rounds a cylindrically coUimated flow. The radius of the 
cylindrically coUimated flow is Rj « Rf/a, where Rf is 
the initial radius of the fast magnetosonic surface. The 
mass flux in the coUimated part is small compared with 
the total mass flux in the wind (~ 10~^), even for rela- 
tively small Lorentz factors. SC essentially decreases for 
larger Lorentz factors. Acceleration of the plasma occurs 
in the radially expanding wind while the plasma in the jet 
is not accelerated in this cold outflow model. 

The relatively small part of the mass flux in the cylin- 
drically coUimated flow has been also found in our simula- 
tions of nonrelativistic outflows (Bogovalov & Tsinganos 
1999 ), although the magnetic flux contained in the rel- 
ativistic case appears smaller by a factor of order 10. In 
addition, the collimation scale at low latitudes does not de- 
pend on the magnetic fleld and the angular velocity of the 
central source for a > 1, in contrast with our results for 
the nonrelativistic winds in the same model. This means 
that the weak collimation of the relativistic plasma is due 
to the intrinsic property of the plasma and not to a prop- 
erty of the model. Collimation decreases due to the affect 
of the electric force which becomes comparable with the 
collimating Lorenz force in the relativistic 

limit and in the high a limit the collimation is satu- 
rated due to the contribution of the electromagnetic fleld 
into the inertial mass of plasma. 

The most striking feature of the coUimated flows ob- 
tained in the model under consideration is that a very 
small fraction (~ 10"'^) of the mass flux exists in the col- 
limated flow in comparison to the mass flux in the wind. 
This result is apparently in conflict with observations of 
AGN. The source of the wind energy is the gravitational 
energy released by the accreted matter. Therefore, the en- 
ergy flux in the wind should be less than the luminosity 
of the accretion disk where about half of the gravitational 
energy is released (Shakura & Sunyaev 1972). The total 
luminosity of AGN Ltot in soft electromagnetic emission ( 
below X-rays) should be comparable with the energy re- 
leased at the accretion. According to our results the energy 
flux in the jet Ejct should be by the factor 10^'^ less than 
the accretion energy, Ejet < 10~^ ■ Ltot- This disagrees 
with observations. X-ray and gamma-ray emission is usu- 
ally a ttributed to the emission from jets (Urry & Padovani 



1995). Sometimes, the luminosity of the jets is less only by 



a factor of 10 than the to tal lu minosity of AGN at other 
wavelengths (Elvis et al. 1994 ). The energy flux in this 
emission can be considered as the lower limit on the en- 
ergy flux in the jet (however, we do not take into account 



why the collimation scale (|36|) at low latitudes does not beaming here). If so, we have in AGN Ejct > 10 -L 



total • 



depend on the magnetic fleld and angular velocity of the 
source. 



This disagreement between theory and observations 
implies that either the present model has no observed 
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counterpart or the SC mechanism is not able in princi- 
ple to provide coUimation of the major fraction of the 
outflow into the jet. It is important in this regard to un- 
derstand from the theoretical point of view what mod- 
ifications should be made in order to provide magnetic 
collimation of the major fraction of the relativistic wind. 
This problem has been already investigated by Bogovalov 



& Tsinganos (2001) in relation to nonrelativistic winds. 



It was shown in this study that in a model consisting of 
a central source ejecting a high-speed wind and an accre- 
tion disk ejecting a nonrelativistic magnetised outflow, all 
the mass flux from the central source can be cylindrically 
collimated indirectly by the surrounding disk- wind. Ap- 
parently, a collimation of the major part of a relativistic 
outflow from a central source can be provided by the same 
mechanism in the context of the present model as well. 
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Appendix A: Details of the numerical simulation 
in the nearest zone. 

For the numerical solution of the time dependent problem, 
the system of equations should be presented as 



dim=200x200 t=1.50 



(A.l) 



where $ is the vector of the plasma state which includes 
the density, components of the velocities and magnetic 
fields. $i and are the first and the second order deriva- 
tives of $ on the spatial variables. Eq. ([|) contains the 
time derivative E in the right hand part generated by the 
displacement current in Eq. (||). Therefore, some special 
transformation of Eq. is necessary to get form (A.l). 
The specific form of this transformation is not unique. In 
our previous work (Bogovalov 1997a| ) unsuccessful choice 
of this transformation resulted in the instability of the nu- 
merical code in the Poynting flux-dominated regime. Here 
we use another form of the equation of motion which pro- 
vides stable numerical code in all regimes. 

This equation can be obtained as follows. Time deriva- 
tive of the frozen-in condition (0) gives that 



E = --(U X B) - -{U X B) - ^E. 

7 7 7 



(A.2) 



Substitution of Eq. ( |A.2| ) and Eq. (y) together with the 
relativistic relationship 7 = ViUi/c in Eq. (|l]) gives the 
equation of motion in the form 

M.kUk = -p(v • V)C/,: + qE,/c + -^((V x B) x B), 

47rc 



1 



where 



K;(B-B)-B,(B-v)), 



(A.3) 



(BiBk + SiVk/c) 



and Sj = (E x B)^. Terms with B in the right hand part of 
Eq. ( [A.3| ) are ehminated with help of Eq. (||). After that, 
multiplication of equation of motion ( [A.3| ) on the inverse 
matrix 



1 



_B2_ 

47rc'^7 



B^Bk , {w ■B/c)S^Bk/A^rpc^-^ + S^Vk/c. 



{Anpc^-f + B^ - (S • v/c) 



(A.4) 




Fig. A.l. The simulation of propagation of the switch on 
shock wave in the relativistic wind. All the notations are 
the same as in Fig. |^ and ^. The wave is generated at the 
beginning of rotation of the central object. 

gives the equation of motion in the form (|A.l| ) appropriate 
for the numerical solution. 

A two-step Lax-Wendroff numerical scheme was used 
for the numerical solution of the problem (Press et al. 
1990|) . The analytical solution obtained for the slow rota- 



tors (Bogovalov 1997a ) was used as the initial state of the 
plasma. 

The testing of the numerical code was performed 
by comparison of the results with the results obtained 
with the nonrclativistic version of the code (Bogovalov 
& Tsinganos 1999| ), by comparison with predictions of the 
analytical theory (Bogovalov 1992, Bogovalov 1997a) and 
by simulation of the perturbation propagation in the wind. 
One of these tests is presented in Fig. A.l. The beginning 
of the star rotation starts the shock wave propagating with 
the Alfvenic velocity in the isotropic wind. This velocity is 
very close to 1 for the taken parameters of the wind. The 
shock has an oscillating structure typical for the collision- 
less shocks. The front is smoothed due to the numerical 
diffusion. The first (largest) maximum can be accepted as 
the boundary of the shock. It follows from the figure that 
the wave propagates with the velocity very close to the 
light velocity which is equal c/vq — 1.01 in the dimension- 
less variables normalized to Vq at 70 = 7. 



